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ABSTRACT 

The solar flare on July 30, 2011 was of a modest X-ray class (M9.3), but 
it made a strong photospheric impact and produced a “sunquake,” observed 
with the Helioseismic and Magnetic Imager (HMI) on NASA’s Solar Dynamics 
Observatory (SDO). In addition to the helioseismic waves (also observed with 
the SDO/AIA instrument), the flare caused a large expanding area of white-light 
emission and was accompanied by substantial restructuring of magnetic fields, 
leading to the rapid formation of a sunspot structure in the flare region. The 
flare produced no significant hard X-ray emission and no coronal mass ejection. 
This indicates that the flare energy release was mostly confined to the lower 
atmosphere. The absence of significant coronal mass ejection rules out magnetic 
rope eruption as a mechanism of helioseismic waves. We discuss the connectivity 
of the flare energy release with the electric currents dynamics and show the 
potential importance of high-speed plasma flows in the lower solar atmosphere 
during the flare energy release. 

Subject headings: Sun: flares; Sun: helioseismology; Sun: sunspots; Sun: X-rays 
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1. INTRODUCTION 


The flare energy release is 
to magnetic reconnection (e.g. 


nsually accornp anied by restrncturing of magnetic field dne 


Su et al.l 


20131) . dissipation, and eruptions of magnetized 


plasma. The restructuring is observed in the form of arcades of magnetic loops with 
reducing shear, changes of the strength and geometry of photospheric magnetic fields, 
erupting magnetic structures (CMEs), and so on. Sometimes, the magnetic topology of 
the flare energy release site can be very complicated due to a complex interaction between 
magnetic fields rooted in multipolar active regions. Despite the complexity, all flares are 


accompanied by standard physical processes like particle accele ration, p 


to several MK and plasma motion. The standard flare model flSturrock 


asma heating up 


19891) assumes 


that the flare energy release processes are caused by magnetic reconnection in the corona, 
but the idea of flare initiation in the lower s olar atmosphere has recently become popular 


because of new observational evidences (e.g. 


Fletcher et ah 


20111) . Physically, the lower 


solar atmosphere significantly differs from the corona. The low atmosphere is characterized 
by partially ionized plasma, higher plasma /9, tu rbulent moti o ns an d highly pronounced 
NLTE conditions during solar flares. The work of[ 


Leake et ah 


(120121 ) is an example showing 


the complexity of magnetic reconnection in the chromospheric plasma and reveals a fast 
reconnection regime with the rate independent on the plasma resistivity. Generally, the 
physical processes in the chromospheric conditions are not well-understood, and will be 
an important topic of future theoretical and observational studies. The current research is 
devoted to the observational study of flare processes in the lower solar atmosphere. 


In this work, as a case study, we present observations and analysis of a moderate M9.3 
class flare, which reveal several significant perturbations in the low atmosphere. The flare 
generated helioseismic waves, caused a large expanding area of continuum emission, and was 
accompanied by substantial restructuring of magnetic fields in the flare region, leading to a 
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rapid formation of a sunspot structure in the flare region. The flare produced no signihcant 
hard X-ray emission, no coronal mass ejection, and no II and III type radio bursts. This 
indicates that the flare energy release was probably conhned in the lower atmosphere in 
closed magnetic structures. The absence of significant coronal mass ejection and II type 
radio burst rules out the magnetic rope eruption as a mechanism of helioseismic waves. 


Helioseismic waves, known as “sunquakes,” can be generated by several mechanisms: 
momentum transf er due to heating of the chrorn ospheric plasma by accele rated energet ic 


charge d particles flKosovichev fc Zharkova 


I 995 II . im pulsive Lorentz force 


Fisher et al. 


(120121) . pressure gradient due t o flux rope eruption fIZharkov et al 


dissipation of electric currents (jSharvkin et al. 


20131) . and rapid 


2014J) . However, t 


sunquakes is still unknown (for a recent review, see 


Kosoviche 


re pr ecise mechanism of 


2014 ). The main goal 


of this work is to trace the plasma dynamics and magnetic held changes in the solar 
atmosphere, as well as possible agents of the hare ene rgy release and sunq uake initiation 


for this particular hare. We use the SDO/HMI data fjScherrer et al. 


20121) to determine 


the hare region magnetic held topology, estimate electric currents, and investigate the 
sunquake . X-ray data fro m the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) 


satellite flLin et al. 


20021) are used to calculate parameters of accelerated particles, and 


estimate their impact in the low chromosphere. 


2. GENERAL CHARACTERISTIC OF THE FLARE 

The M9.3 hare on July 30, 2011, occurred in active region NOAA 11261, which had a 
J-type magnetic conhguration. According to the GOES-15 soft X-ray data (Fig. 1), the 
hare started at 02:04 UT, reached maximum at 02:09 UT, and ended at about 02:30 UT. 
The RHESSI observations (Fig. 1) show the presence of the X-ray emission up to 300 keV, 
reaching maximum at 02:08 UT. 
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According to the Learmonth and Cnlgoora radio observatories, there were no type II 
radio bnrsts. This means the absence of signihcant shock waves, which conid be associated 
with ernpting magnetic strnctnres. There was no recorded CME that indicates on a 
non-ernptive natnre to the event. The radio data also show the absence of III type radio 
bnrsts, thns indicating a closed geometry of magnetic held lines in the haring region. 


The Atmospheric Imaging Assembly (AIA) on SDO flLemen et al 


2 OI 2 I) observes the 


Snn with time cadence 45s spatial resolntion 1.2" (pixel size 0.6"). In Fig. 2 we show 
the AIA images before and after the hare, illnstrating the reconhgnration of ribbon-like 
strnctnres in the AIA 304 A channel {IgT = 4.92). This reconhgnration probably 
corresponds to a redaction of the magnetic shear dnring the hare. In the AIA 94 channel 
A {IgT = 6.86), we observe the formation of hot large-scale loops. These AIA images are 
taken at the prehare and posthare times. 


The hare signals located in the vicinity of the magnetic held polarity inversion 
line (PIL) are also detected in all HMI observables (Fig. 3). We will discnss the HMI 
observations in the next sections. 


3. CONTINUUM EMISSION AND FLOW PATTERN IN THE FLARE 

REGION 

The hare prodnced a strong emission observed in the HMI data continnnm data. It 
originated near PIL and expanded in the North-West direction in a wave-like manner 
(Fig. 3b). The peak of the HMI continnnm intensity averaged over the hare region has a 
2-minnte delay relative to the HXR maximnm (Fig. 1-top). This continnnm signal is shown 
in Fig. 3b by the time-diherence images. The estimated velocity of the signal expansion is 
~ 30 — 40 km/s, which is higher than the local sonnd speed and can be comparable with 
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the magneto-acoustic speed of ~ 20 km/s. 

We also show in Figs. 3b and 3c the line-of-sight velocity structure by contours. 

The velocity variations are revealed in the Dopplergram differences (Fig. 3a). The most 
intriguing feature is the strong upflows reaching 3.6 km/s in the vicinity of PIT. These flows 
are long-lived structures seen before the flare, so they cannot be artifacts due to impulsive 
variations of the line prohle caused by the flare. Besides, we observe sharp changes in the 
flow structure during the most intensive HXR emission. Such impulsive variations may be 
associated with the rapid variations of the line profile, and in this case, the corresponding 
Doppler velocity variations are difficult to estimate. After the flare impulsive phase, the 
upflow region elongated along PIT is divided into two upflow regions separated by a distance 
of ~ 4 Mm. 


4. MAGNETIC RESTRUCTURING 

The flare is accompanied by significant restructuring of the magnetic field, which 
resulted in a rapid formation of a small sunspot from an initially diffuse magnetic field. The 
spot formation happened during the short flare period < 1 hour (Fig. 4). To our knowledge, 
this is the hrst observation of such rapid sunspot formation associated with a solar flare. 

Vector magnetograms from HMI also reveal substantial changes of the magnetic field 
structure in the vicinity of the magnetic neutral line during the flare. In particular, the 
magnetic shear is sharply reduced (Fig. 5). To illustrate the magnetic field reconfiguration. 


we calculate a non-linear force-free field extrapolation (NLFFF 

) of magnetic field using the 

) with 720 s time cadence 

disambiguated HMI vector magnetograms ( 

Centeno et ah 

2014 

as the initial condition. We use the technique described in ('' 

Vheatland et ah 

o 

o 

o 

C^l 


based on minimization of functional Jy\j x B\‘^ + |V ■ B\‘^dV. The results of the NLFFF 
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extrapolation for the vector magnetograms before and after the flare are shown in Fig. 5. 
Initially, a sheared magnetic structure elongated along the PIL (marked by red color) 
became less sheared after the flare. This corresponds with a transition to a more potential 
magnetic structure. It appears that these changes in the magnetic topology caused a rapid 
concentration of the magnetic held and formation of the small sunspot (Fig. 4). 


Previous observa tions show ec 
after some hares (e.g. IWang et ah 


the disappearance of small sunspot-like structures 


2002 


Chen et ah 


20071) . but this event is probably the 


hrst observation of the formation of such a magnetic structure caused by hare magnetic 
restructuring. 


5. SUNQUAKE 


The sunquake event was initially revealed as a circular shaped wave in the running 
diherence of the Doppler velocity data. The sunquake wave (Fig. 6) is best visible about 20 
minutes after the initial hare impact on the photosphere located near the PIL (Fig. 3). The 
wave front was highly anisotropic. The wave front traveling outside the magnetic region to 
the South-East had the highest amplitude. In the North-West direction, the wave traveled 
towards a big sunspot, and its amplitude was suppressed when it reached the sunspot. 


Figure 6b illustrates the positions of the hare impact and the helioseismic fronts 


observed in the Dopple r-shift data, hltered for a cen 


time-distance diagram flKosovichev fc Zharkova 


1998 


ral frequency of 6 mHz. The 


Zharkova fc Zharkov 


29971) for the 


helioseismic wave traveling toward the South-East is shown in Fig. 6a. The dashed curve 
represents the theoretical time-distance relation calculated in the ray approximation for 
a standard solar model. The location of this curve is chosen to approximately match the 
leading wave front. The short strong signal at the beginning of the event indicates the 





















initial source motion. 


The time-distance diagrams show two interesting features. First, during approximately 
the hrst three minutes, the wave source moved with a speed of about 15-17 km/s, which is 
higher than the local sound speed but may correspond to the magneto-acoustic speed in the 
sunspot penumbra in the vicinity of the source. Second, this motion may be associated with 
a supersonic plasma eruption, or caused by a series of impacts due to a fast moving magnetic 
reconnection process. The source motion is a common feature of sunqu akes, and may be 


essential for generating high-amplitude helioseismic waves in solar flares flKosovichev 


201 4h 


The strongest waves are usually observed in the direction of the source motion flKosovichev 


200d). 


6. ELECTRIC CURRENTS IN THE FLARE REGION 

In this section, we consider the evolution of the electric currents at the photosphere 
level. To estimate the horizontal electric currents, we use Faraday’s law applied to the 
45-second line-of-sight HMI magnetograms with the spatial resolution V and pixel size 0.5": 


E-dl 


c 


c dt 



( 1 ) 


We can estimate the average transversal component of electric held {E±) = [d^z/dt]/cL, 
where is the total magnetic hux inside a contour with length L, which covers the hare 
region. The evolution of d^z/dt presented in Fig.l (histogram in bottom panel) shows that 
the whole hare impulse correlates with d^z/dt. 


To calculate the verti cal currents we use the d isambiguated H 


data flCenteno et ah 


201411 . and the Ampere’s law (iGuo et al.l 


Ml vector magnetic held 


2018 ll 
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3 z 


Uv X B), 


c_ 

47r \ dy dx ) 


( 2 ) 


The resulted maps, effectively averaged over 12 min due to the HMI temporal 
resolution, are presented in Fig. 7. We can see that the enhancement of electric currents 
partially correlates with the location of the strongest photospheric impact. 

To estimate the measurement error, we calculate the distribution of electric currents in 
a non-flaring region, and then approximate by a Gaussian. The resulted standard deviation 
is considered as error of jz due to noise. The time evolution of (j^) averaged over the 
flare region (Fig. 8) shows that the maximum of {jz) is significantly delayed (~ 30 min) 
relative to the X-ray flux observed from GOES. Perhaps, the generation and dissipation of 
electric currents continued after the flare impulsive phase, and might be associated with the 
magnetic restructuring. It might also be explained by not accounting for the fine structure 
of electric currents in the HMI data. The problem is we assume that the observed current 
density is distributed uniformly through HMI pixels, but it could be concentrated in much 
thinner tubes. Assuming small values of the filling factor at the flare beginning and its 
subsequent increase due to the current dissipation, one could fit (j^) to the X-ray time 
profile. Resolving this issue requires high-resolution spectro-polarimetric observations. 


7. RHESSI IMAGES AND SPECTRA 


To determine properties of the accelerated particles and the hot flare plasma, we 
use the RHESSI data in the range of 5-250 keV. RHES SI employs a Fourier technique to 


reconstruct X-ray emission sources 


Hnrford et al 


( 20021 ). We applied the GLEAN algorithm 


to synthesize the X-ray images. In Fig. 9 (right panel), the RHESSI HXR and SXR 
contour images are compared with the corresponding AIA 94 A images for the time interval 
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covering the HXR peaks of both flares. We choose this AIA channel because it experiences 
less saturation during the impulsive phase. The AIA brightnings correlate with the HXR 
emission sources and reveal a loop-like structure with a more intensive West footpoint. 


The power-law approximation f{E) = AE (where A is a normalization coefficient 
for energy 10 keV) is considered for the hard X-ray (HXR) nonthermal emission 20 keV . 


To si mulate the presence of the low-energy cutoff, we use the broken power law flHolman 


20031) with hxed photon spectral index 70 = 1.5 below the break energy {Eiow] 


The thermal soft X-ray (SXR) spectrum at < 20 keV is approximated by a single¬ 
temperature thermal bremsstrahlung emission law with two parameters: temperature (T) 
and emission measure {EM). The RHESSI spectra are htted by means of the least squares 
technique, implemented in the OSPEX package with 5 free parameters {EM, T, A, Eio^,, 
and 7 ). Figure 9a (top) displays the observed X-ray spectra and the results of the model 
htting with = 2.37. 


Th e total flux of accelerated elect rons, E [electrons s 


-11 


work of 


Svrovatskii fc Shmeleva! (19721): 


can be estimated following the 


E{Eiouj < E < Ehigh) = 1.02 X 10 


34 


Iph{EloW ^ E <i Ehigh) 


(3) 


Eiow/3{Si, 1 / 2 ) [1 — {Eiow/Ehigh)^^] 

where Ehigh is the upper energy cutoff, 5 = 7 — 1 is the spectral index of accelerated electrons 
in the HXR emission region, (3{x,y) is the beta function, and Iph{Eio^ < E < Ehigh) 
photons s“^ cm“^ is the energy of the integrated photon spectrum in the range shown 
in the brackets. From the htting results using this formula, we obtained electron hux 
E ^ (1.6 ± 0.4) X 10^^ electrons s“^. The total energy hux of accelerated electrons is 
Pnonth = 4.6 X 10^® erg s“^. Theoretically, these electrons could contribute to the sunquake 


initiation. However, the discrepancy between the locations of the sunquake impact and the 
strongest HXR emission source indicates that the beam-driven origin of the sunquake is 








unlikely. 


8. DISCUSSION 


The current paradigm of solar flares is that the magnetic energy is released in 
the solar corona and transported towards the lower atmosphere by heat and energetic 


particles. For instance, in the hvdrodvnamic thick-tareet model fe.e. 

Kostiuk & P 

ffilner 

1975; 

Fisher et al.l 

1985; 

Kosovichev 

1986; 

Allred et al. 

2006; 

Rubio da Costa et al. 

2014). 


high-energy electrons accelerated in the upper corona are injected along magnetic held 
lines into the atmosphere, generate hard X-ray emission in the loop footpoints, and heat 
the upper chromosphere to high temperature, producing a high-pressure region. The 
high-pressure region expands producing upward and downward propagating shocks. The 
downward shock may reach the photosphere and cause a sunquake. 


However, this M-class hare showed a very strong photospheric and helioseismic signal, 
but had relatively weak high-energy emissions and coronal dynamics. These observations 
suggest that a substantial part of the magnetic energy is released directly in the low 
atmosphere. This requires a new mechanism of energy release and transport into the low 
atmosphere during the hare impulsive phase. One of the possible alternative mechanisms 
of sunquake and hare initiatio n can be a rapid dissipation of electric currents or sharp 


enhancement of Lorentz force f|Fisher et al 


2012 


Sharvkin et al. 


20141) . The dynamics and 


spatial structure of electric currents give evidence in favor of such a hypothesis. However, 
the relationship between the current dynamics and precesses of particle acceleration is 
unclear. The observations show that the maximum X-ray huxes occurred before the mean 
electric current reached maximum. 


The observed wave-like continuum emission probably is associated with the sunquake 
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initiation. The direction of propagation of this wave corresponds to the wave observed in 
the HMI LOS velocity. The estimated velocity is ~ 30 — 40 km/s, which is higher than 
photospheric magneto-acoustic speed \/v\ + ~ 20 km/s for a 500 G magnetic held. 


The observed uphows across the magnetic held near the PIL may play a very important 
role in hare initiation process. These hows could lead to a transversal electric held 
Et = \v X B\/c ~ VzBt/c, generating circulating electric currents between the solar lower 
atmosphere and the corona (Fig. 10). In the conditions of partially ionized magnetized 
plasma, ambipolar dihusion can lead to efficient dissipation of impulsive electric currents, 
since coefficient may be seve ral orders higher than the coefficient of Ohmic dihusion (e.g. 


Khomenko fc Colladoi 


20121) . To clarify such an idea, a detailed modeling of the hare 


plasma dynamics and heating using a multi-huid MHD approach is needed. One of the main 
peculiarity of the considered hare is a fast formation of the small sunspot. The magnetic 
held evolution is described by induction equation: 


dB 


V X {v X B) + fiAB 


where fi = c^/(47rcr) is magnetic dihusivity for plasma with conductivity a. In the 
case of our hare observed sunspot formation can not be connected with dihusion term as 
it will lead only to sunspot breakup. Probably hows is the main agent of such magnetic 
held restructuring. Observed by the HMI vertical uphows probably will contribute to the 
magnetic held intensihcation as dx{vzBx) + dy{vzBy). However, in the case of turbulent 
medium v = vq + 5v and B = Bq + 5B with averages |5n| = \5B\ = 0, time averaged 
induction equation has form: 


dBo 


dt 


V X (uq X Hq) “1- V X {5v X SB'j -|- 
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where ^turb is turbulent magnetic diffusivity which is larger than classical one. Thus, 
we have three competing processes of magnetic held changing in the hare region: convective 
generation of magnetic held by laminar (1st term) or stochastic (2nd term) hows and 


magnetic dihusion (3d term). 


flWang et al 


2002 


Chen et al. 


n a 


ot of hares we usually observe breakup of sunspots 


20071 ). probably connected with fast dihusion process 


(magnetic reconnection), and thus, the convective terms are not very important. In our 
event the hows generating magnetic held can be more important than dihusion and lead 
to the formation of a small sunspot from initially dihuse sunspot region. To conhrm such 
ideas and observations we need detailed self-consistent numerical MHD simulations. 


9. CONCLUSIONS 

The main results of the work are the following: 

1. The SDO/HMI observations reveal signihcant changes of the magnetic held topology, 
leading to a rapid sunspot formation during the solar hare. 

2. The HMI continuum emission with velocity of ~ 30-40 km/s shows a fast wave-like 
expansion. This wave may be related to a shock generated by the energy release in 
the low atmosphere. 

3. The HMI Doppler shift data detect a persistent plasma uphow across the magnetic 
held near the polarity inversion line with a subsequent break into two uphow regions 
after the hare impulsive phase. The uphow may play an important role in maintaining 
the circulating electric currents contributing to the energy release. 

4. The location of the strongest photospheric disturbance leading to the sunquake 
does not correlate with the HXR emission source. This is not consistent with the 
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hypothesis of the standard flare model. The sunquake generation place is located in 
the vicinity of the npflow break up, and partially corresponds to the most intense 
electric currents. 

The presented observations show that the flare energy release in the low solar 
atmosphere not only associated with the precipitation of high-energy particles from the 
corona. The detected electric currents may also be responsible for generating the flare 
event. However, the understanding of the mechanism which causes rapid current dissipation 
requires high-resolution observations and theoretical modeling. 

The work was partially supported by RFBR grant 13-02-91165, President’s grant 
MK-3931.2013.2, NASA grant NNX14AB70G, and NJIT grant. 
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Fig. 1.— Top: color curves show the RHESSI count rates in four different energy bands 
(marked in the hgure); the histogram shows the HMI averaged continuum intensity. Bottom: 
GOES light curves in two channels 0.5-4 A (red) and 1-8 A (black); the histogram shows the 
time derivative of the magnetic flux calculated from the HMI data. 
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Fig. 2.— The SDO/AIA images in two channels (304 A, top) and (94 A, bottom), before 
(left) and after (right) the flare. Red contour curves show the polarity inversion line. 
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Fig. 3.— HMI images: a) time differences of Dopplergrams, b) time differences of HMI 
continnnm intensity images, and c) seqnence of HMI continnnm intensity images. Black 
line shows the polarity inversion line. Red and bine contonrs show regions of npflows and 
downflows according to the HMI Doppler measnrements. 
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Fig. 4.— Sequence of the HMI continuum intensity images. Black line shows the polarity 
inversion line. Red and blue contours show regions of upflows and downflows according to 
HMI Doppler measurements. 



Fig. 5.— NLFFF reconstruction of magnetic held. Left and right panels correspond to 
prehare and posthare times. Red lines are marked to show the evolution of low-lying magnetic 
held lines near the magnetic polarity inversion life. 
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Fig. 6.— (Left) The sunquake time-distance diagram with the ray-path theoretical predic¬ 
tion (dashed yellow line); (Right) locations of the sunquake source and waves (marked by 
yellow arrows). 
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,— Sequence of 12 images which shows the evolution of the vertical component of 
current density (color bar from -0.13 to 0.17 A/m^). Black line shows the polarity 
)n line; white and gray contours show upflows and downflows according to the HMI 
r measurements. Red contours show the LOS magnetic held at 1,000 and 1,500 G 
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Fig. 8.— Time evolution of the mean vertical current density (red curve shows (j^) > 0, 
black curve shows (jz) < 0). Blue line is the GOES 1-8 A lightcurve. 
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Fig. 9.— a): RHESSI X-ray spectra (upper panel), with the htting results printed in the 
plot, the htting residuals are shown in the bottom panel, b): AIA 94 A image with red and 
blue contours showing the X-ray sources with energies 6-12 and 50-100 keV. 
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Fig. 10.— The scheme of the flare region near the polarity inversion line. 








